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Abstract 

We study the synchronization of a coupled map lattice consisting of a one-dimen- 
sional chain of logistic maps. We consider global coupling with a time-delay that 
takes into account the finite velocity of propagation of interactions. We recently 
showed that clustering occurs for weak coupling, while for strong coupling the array 
synchronizes into a global state where each element sees all other elements in its 
current, present state [Physica A 325 (2003) 186, Phys. Rev. E 67 (2003) 056219]. 
In this paper we study the effects of in-homogeneities, both in the individual maps, 
which are non-identical maps evolving in period-2 orbits, and in the connection links, 
which have non-uniform strengths. We find that the global synchronization regime 
occurring for strong coupling is robust to heterogeneities: for strong enough average 
coupling the inhomogeneous array still synchronizes in a global state in which each 
element sees the other elements in positions close to its current state. However, the 
clustering behaviour occurring for small coupling is sensitive to inhomogeneities and 
differs from that occurring in the homogeneous array. 
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Globally coupled units with time delays and inhomogeneous interactions arise 
in a variety of fields. In the fields of economics and social sciences, a multi-agent 
model that allows for temporally distributed asymmetric interactions between 
agents was recently proposed in Ref. [1]. The model defined a coupled map 
lattice with interactions obeying a Gaussian law and transmitted through 
a gamma pattern of delays. In the field of optics, a globally coupled laser 
array with feedback has been proposed for exploring the complex dynamics of 
systems with high connectivity [2,3,4]. Globally and locally coupled maps have 
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been employed to study synchronization [5] in a great variety of fields ranging 
from activity patterns in pulse-coupled neuron ensembles [6,7] to dynamics 
models of atmospheric circulation [8]. 

Since they were introduced by Kaneko, globally coupled logistic maps have 
turned out to be a paradigmatic example in the study of spatiotemporal dy- 
namics in extended chaotic systems [9]. In the simplest form all maps are 
identical and interact via their mean field with a common coupling: 

e ^ 

x.{t + !) = (!- e)f[x.{t)] + -Y: f[xM- (1) 

j=i 

Here t is a discrete time index, i is a discrete spatial index: i — 1 . . .N where 
N is the system size, f{x) = ax{l — x) is the logistic map, and e is the coupling 
strength. Even though the model has only two parameters (the common non- 
linearity a and the coupling strength e), it exhibits a rich variety of behaviours. 
For a large e the maps synchronize in a global state, and evolve together as 
a single logistic map. For intermediate coupling the array divides into two 
clusters {xi{t) = x.j{t) if i and j belong to the same cluster) which oscillate 
in opposite phases. For small coupling the number of clusters increases but 
is nearly independent of the total number of maps in the system. Finally, for 
very small coupling the number of clusters is proportional to N [10]. 

We have recently studied effect of time-delayed interactions [11,12]. We con- 
sidered the array: 

e ^ 

x,{t + 1) = (1 - e)f[x,{t)] + -Y: f[xAt - r^J)] (2) 

where Tij is a time delay, proportional to the distance between the ith and 
jth maps. We took Tij = k\i — j\, where k is the inverse of the velocity of the 
interaction signal traveling along the array. 

We found that for weak coupling the array divides into clusters, and the be- 
havior of the individual maps within each cluster depends on the delay times. 
For strong enough coupling, the array synchronizes into a single cluster (glob- 
ally synchronized state). In this state the elements of the array evolve along 
a orbit of the uncoupled map, while the spatial correlation along the array is 
such that an individual map sees all other maps in his present, current, state: 

Xj{t~Tij) =Xi{t). (3) 

It was also found that for values of the nonlinear parameter a such that the un- 
coupled maps are chaotic, time-delayed mutual coupling suppress the chaotic 
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behavior by stabilizing a periodic orbit which is unstable for the uncoupled 
maps. 



In this paper we extend the previous study to assess the influence of hetero- 
geneities. We consider the following array: 

1 ^ 

Xi{t + 1) = (1 - ei)f[ai, Xi{t)] + ^Y1 ^ij/K-' ^A't - nj)] (4) 



where f{a,x) = ax{l — x), Sij is the strength of the link coupling the maps i 
and J, Si is the average connection strength of the site i: 

1 ^ 



and Tij = k\i—j\. We perform numerical simulations with the aim of studying 
the effects of inhomogeneities in the maps and in the strength of the couphng 
hnks. We consider two different situations: 

1) Nonidentical maps coupled with identical coupling strengths (ejj = e\/ i, j). 
The value of a, is random, uniformly distributed in the interval {ao—5a,ao+Sa). 
We limit ourselves to consider maps in the period-2 region (oj G [3, 3.449 ...]). 
The study of the synchronization of the array when the individual maps with- 
out coupling are either in fixed points or in period-n orbits is left for future 
work. 

2) Identical maps (oj = V i) coupled with nonindentical links. In this case 
we consider two different situations: (a) The value of Sij random, uniformly 
distributed in the interval (0,e). (b) Sij varies in interval [0,e], such that the 
strength of a link decreases with its length. We take Sij = e[l — Tij/ max(rij)]. 
It can be noticed from Eqs. (4,5) that the globally synchronized state Eq.(3) 
is also an exact solution of the array when the maps are identical but the 
connection strengths are not. 

Let us first illustrate the transition to global synchronization in the case of 
inhomogeneous maps connected through homogeneous links (ej.j = e V i. j). We 
exemplify results for k odd. For k odd the globally synchronized state of the 
homogeneous array is the anti-phase state: xf{t) = Xa, Xf,, Xa, ■ ■ ■; xf{t -|- 1) = 
xi,,Xa,Xb, . . .. Here Xa and Xb are the points of the limit cycle of the map 
f{ao,x). Similar results are observed for k even (in this case the globally 
synchronized state is the in-phase state: xj{t) — Xa,Xa,Xa, ■ ■ ■', x^{t -|- 1) = 

Xh, Xh, Xf), . . .) . 

Figures 1(a) and 1(b) display bifurcation diagrams for increasing e, which 



3 




Fig. 1. Bifurcation diagram for homogeneous (a) {6a = 0.) and inhomogeneous (b) 
{6a = 0.2) array. Parameter values: a = 3.2, k = 1, and N = 100. 
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Fig. 2. Array configuration at time t = 3000, for different values of the coupling 
strength e (a) e = 0.07, (b) e = 0.11, (c) e = 0.29, and (d) e = 0.98 . Other 
parameters as in Fig. 1 . 

were done in the following way: we chose a random initial condition [xi{0) 
with i — 1, . . . ,N], and plotted a certain number of time-consecutive values of 
Xi{t) (with t large enough and i = 1, . . . , A^) vs. the coupling strength e. Since 
the initial condition fixes the cluster partition of the array, the same initial 
condition was used for all values of e. For comparison. Fig. 1(a) displays the 
bifurcation diagram for a homogeneous chain {5a = 0), while Fig. 1(b) displays 
results for 6a — 0.2. While a sharp transition to global synchronization is 
observed for 5a — 0, a, smoother transition, reminiscent of bifurcations in the 
presence of noise, is observed for 6a ^ 0. To illustrate this transition with 
more detail we show in Fig. 2 the configuration of the array for several values 
of e (the values of e are indicated with an arrow in Fig. 1). It is observed that 
the array gradually evolves to the antiphase state. For low e there is a large 
dispersion in the values of xi and defects are observed. The dispersion and the 
number of defects gradually diminishes as e increases. 

Next we study the case of homogeneous maps (oj = Oq V i) connected through 
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Fig. 3. Bifurcation diagram displaying the transition to the globally synchronized 
state. N = 100, Oj = 3.2 V i, and k = 1. (a) e^j = e^y i and j; (b) e^j is uniformly 
distributed in (0, e) ; (c) £j j decreases with distance. 

inhomogeneous links. Figure 3 displays bifurcation diagrams done in the same 
way as in Fig. 1, i.e., taking the same initial configuration of the array for all 
values of e. However, in this case the links have different coupling strengths 
[Sij varies in (0,e), is either uniformly distributed or decreasing with Tij]. 
Therefore, in the horizontal axis we plotted the average coupling strength 
< e >— l/N'^J2iJ2j£i,j- If can be observed that as < £ > increases there 
is a sharp transition to global synchronization, both when eij is randomly 
distributed in (0,e), and when e^j is in (0,e) decreasing with distance. However, 
there is different clustering behaviour, as illustrated in Fig. 4. Figures 4(a) and 
(b) display the clustering behaviour for /c = 1 (we point out that for /c = 1 
and < £ > large enough the array synchronizes in anti-phase), and Figs. 4(c) 
and (d) display the clustering behaviour for A; = 2 (for k — 2 and < e > large 
enough the array synchronizes in-phase). For comparison. Figs. 4(a) and 4(c) 
display the clustering behaviour of the homogeneous array. In these figures 
%j = Co with eo equal to the average coupling strength in Figs. 4(b) and 
4(d). In spite of the fact that the average connection strength is the same, the 
partition of the array is different when the connection links have non-uniform 
strengths. 

To quantify the effect of in-homogeneities we introduce the following quantity 



a(t) 



N 



Ek.(t)-a;f(t)]2, 



(6) 



that measures the deviation of the array configuration, Xi{t), from a given 
reference state, x^{t). Xi{t) is evaluated at time t with t long enough to 
let transients die away. The reference state is the anti-phase state for k odd 
(xf (t) = Xa, Xb, Xa, ■ ■ .; xf{t -|- 1) = x^, Xa, X),, . . .) and the in-phase state for 
k even {x^{t) — Xa, Xa, Xa, ■ ■ .; x^{t -|- 1) = Xb, Xf,, Xf,, . . .). Here Xa and Xf, the 
points of the limit cycle of the map /(oq, x). 

Figure 5 shows < cr > vs. e where < • > represents an average over different 
initial configurations Xi{0). We observe here that in the homogeneous array 



5 




'O 20 40 . 60 80 100 '^'^o 20 40 60 80 100 

1 I 



Fig. 4. Clustering behaviour for weak coupling strength. The thin line shows the 
array configuration at time t, and the thick line shows the array configuration at 
time t + 1. iV = 100 and a,; = 3.2 V i. (a) fc = 1 and e^j = 0.08 V {h) k = 1 and 
Ejj uniformly distributed in [0,0.16], < £ij >= 0.08; (c) k = 2 and Sij = 0.073 V 
(d) k = 2 and '3. 
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Fig. 5. Deviation of the array configuration vs coupling strength averaged over 100 
realizations for the cases with (continuous line, 6a = 0.2) and without dispersion 
(dashed line, Sa = 0). (Parameter values: N = 100, k = 1, a = 3.2) 

(solid line) the decay of the deviation a is more abrupt that in the hetero- 
geneous array (dashed line). This is due to related to the sharp transition to 
global synchronization shown in Fig. 1. 

To summarize, we have studied the effect of the inhomogeneities in a linear 
chain of globally coupled logistic maps with time-delayled interactions. The 
inhomogeneities occurs both in the individual maps which are non-identical 
and in the connection link which have non-uniform strengths. We considered 
the case in which the maps, without coupling, evolve in a limit cycle of pe- 
riod P — 2. We found that the global synchronization regime occurring for 
strong coupling is robust to heterogeneities: for strong enough average cou- 
pling the inhomogeneous array still synchronizes in a global state in which 
each element sees the other elements in positions close to its current state. 
However, the clustering behaviour occurring for small coupling is sensitive to 
inhomogeneities and differs from that occurring in the homogeneous array. 
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